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jlLi' We analyze the accelerator constraints on the parameter space of the Minimal Supersym- 

r-| , metric extension of the Standard Model, comparing those now available from LEP II and 
anticipating the likely sensitivity of Tevatron Run II. The most important limits are those 
^ ■ from searches for charginos x^i neutralinos Xi ^i^nd Higgs bosons at LEP, and searches for stop 
a I squarks, charginos and neutralinos at the Tevatron Collider. We also incorporate the con- 
straints derived from 6 — S7 decay, and discuss the relevance of charge- and colour-breaking 
minima in the effective potential. We combine and compare the different constraints on the 
Higgs-mixing parameter /i, the gaugino-mass parameter mi/2 and the scalar-mass parameter 
mo, incorporating radiative corrections to the physical particle masses. We focus on the 
resulting limitations on supersymmetric dark matter, assumed to be the lightest neutralino 
X, incorporating coannihilation effects in the calculation of the relic abundance. We find 
that rriy^ > 51 GeV and tan/? > 2.2 if all soft supersymmetry-breaking scalar masses are uni- 
versal, including those of the Higgs bosons, and that these limits weaken to > 46 GeV 
and tan/5 > 1.9 if non- universal scalar masses are allowed. Light neutralino dark matter 
cannot be primarily Higgsino in composition. 



1 Introduction 



The search for experimental evidence for supersymmetry is currently approaching a tran- 
sition. For several years now, many of the most incisive experimental searches have been 
those at LEP whose constraints on the parameter space of the minimal supersymmetric 
extension of the Standard Model (MSSM) have grown ever more restrictive, as the centre-of- 
mass energy of LEP II has been increased in successive steps. In parallel, improved analyses 
of data from Run I of the Tevatron Collider have been providing important complementary 
constraints [@|. The transition is marked by the termination of the LEP II experimental 
programme in late 2000 and the anticipated start of Run II of the Tevatron Collider in 2001. 

The results of experimental searches for different MSSM particles can usefully be com- 
pared and combined using the conventional parameterization of the model in terms of 
supersymmetry-breaking scalar and gaugino masses rriQ, mi/2, the higgsino mixing parame- 
ter /X, the ratio of Higgs vacuum expectation values (vev's) tan f3 and a universal trilinear 
supersymmetry-breaking parameter A. We work in the framework of gravity-mediated mod- 
els of supersymmetry breaking, in which it is commonly assumed that the scalar masses 
mo and the gaugino masses mi/2 are universal at some supersymmetric GUT scale. The 
assumptions that these supersymmetry-breaking parameters are universal should be ques- 
tioned, particularly for scalar masses and especially those of the Higgs supermultiplets, but 
provide a convenient way of benchmarking comparisons and combinations of different ex- 
perimental searches. In this paper, we make such comparisons and combinations in variants 
of the MSSM in which the scalar-mass universality assumption is extended to Higgs fields 
(UHM, also commonly referred to as mSUGRA or the constrained MSSM or CMSSM), and 
also without this supplementary assumption (nUHM). 

In making such comparisons, we emphasize the importance of including radiative cor- 
rections to the relations between these MSSM model parameters (mo, mi/2, yU, tan/5, A) and 
the physical masses of MSSM particles. Radiative corrections are well-known to be crucial 
in the MSSM Higgs sector, but also should not be neglected in the chargino, neutralino, 
gluino and squark sectors. As we have emphasized previously 0, the differences between 
the domains of MSSM parameter space apparently explored at the tree and one-loop levels 
are comparable to the differences between the domains explored in successive years of LEP 
running at higher centre-of-mass energies. In view of the intense experimental effort put into 
sparticle searches at LEP II, it is important that the final results of these efforts be treated 
with the theoretical care they deserve. This issue is also relevant if one wishes to compare the 
physics reaches for electroweakly-interacting sparticles at LEP and for strongly-interacting 
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sparticles at the Tevatron Collider, in which case one should take into account the important 
radiative corrections to squark and gluino masses 0], as well as to their production cross 
sections. 

In addition to direct searches for the production of MSSM particles, important indirect 
constraints must also be taken into account. These include other accelerator constraints, such 
as the measured value of the 6 — 57 decay rate |H, ^, and non-accelerator constraints related 
to the possible role of the lightest supersymmetric particle (LSP) as cold dark matter (CDM). 
The lightest supersymmetric particle (LSP) would be stable in any variant of the MSSM 
which conserves R parity, as we assume here. In gravity-mediated models of supersymmetry 
breaking, the framework adopted here, the LSP is commonly thought to be the lightest 
neutralino and calculations of the cosmological relic density of LSPs, Q^, yield values in 
the range preferred by cosmology in generic domains of MSSM parameter space 0. The 
possibility of supersymmetric CDM provides one of our principal motivations for seeking a 
deeper understanding of the allowed MSSM parameter space, but is not our only focus in 
this paper. 

The most essential dark-matter constraint is that the relic LSP density not overdose the 
Universe. The conditions that the universe has an age in excess of 12 billion years and that 
^totai < 1 imply an upper bound on ^^^"^ O-^- Further, the convergent indications from 
astrophysical structure-formation arguments and observations of high-redshift supernovae 
are that ^cdm < 0.5|^, whereas the Hubble expansion rate Hq = lOO/i km/s/Mpc: h = 0.7 
with an error of about 1O%0, so we require VLisph^ < ^CDWih'^ < 0.3. On the other 
hand, astrophysical structure formation seems to require flcDM > 0.2, so we also require 
^Lsph"^ > 0.1, while acknowledging that a lower value of Qlsp could be permitted if other 
CDM particles such as axions and/or superheavy relics are present. 

There have recently been some significant developments in the analysis of supersymmet- 
ric CDM. One is that the importance of co-annihilation effects involving next-to-lightest 
supersymmetric particles (NLSPs) such as the f, fi and e for calculations of the relic density 



of a gaugino-like LSP has recently been recognized |]10|, |TT|. Another phenomenon whose 



importance in the CMSSM has recently been underlined is the possible transition of the 
electroweak vacuum into a charge- and colour-breaking (CCB) minimum [jl2[. The absence 
of such an instability is not absolutely necessary, since a transition in the future cannot be 
excluded. Therefore, we comment on the regions of MSSM parameter space in which the 
CCB instability is absent, but do not focus exclusively on these regions. 

The main purpose of this paper is to prepare for the compilation, comparison and com- 
bination of the definitive results from LEP II and the Tevatron. We illustrate our analysis 
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with the latest available limits from these two experimental programmes |1|, 0], supplemented 
by educated guesses at their final sensitivities. As we have explained previously, and discuss 
in more detail below, a key role in constraining the MSSM parameter space is provided 
by the LEP Higgs search. We express our results as a function of the present LEP lower 
limit on mn, currently 107.9 GeV [|r^, and the prospective future sensitivity, which may 
approach 112 GeV. We use our analysis to present lower limits on the LSP mass and on 
tan/3. We include a discussion of the implications of relaxing the UHM assumption that the 
soft supersymmetry-breaking contributions to Higgs masses are also universal. In particular, 
we investigate whether a light Higgsino LSP is still a viable dark matter candidate, and 
find that the latest LEP II data now exclude this possibility. Finally, we discuss the likely 
future developments in the exploration of the MSSM parameter space in the period before 
the start-up of the LHC, during which the central role is likely to be played by Run II of the 
Tevatron Collider. 

The layout of this paper is as follows. In Section 2 we review in more detail the theo- 
retical framework we adopt, discussing the issues of universality and relic coannihilations, 
and stressing the importance of Higgs mass constraints. In Section 3 we review our imple- 
mentation of the 6 — s> S7 constraint, including, where applicable, the next-to- leading-order 
(NLO) QCD corrections. We discuss the implications of the latest available constraints from 
LEP II in Section 4, combining them in Section 5 with the cosmological and astrophysical 
constraints 0.1 < QcDAih"^ < 0.3 as well as the b ^ s'j constraint, and making the UHM 
assumption. We find 

> 51 GeV, tan/5 > 2.2 (1) 

and discuss the expanded ranges of and tan/5 that may be explored by the improved 
Higgs-mass limits that might be obtained from the run of LEP II in the year 2000. The 
limits in (|T]) are strengthened when we restrict values of Aq to minimize the parameter space 
with CCB minima, in which case we find 

> 54 GeV, tan/5 > 2.8 (2) 

We further generalize the discussion to non-universal Higgs masses (nUHM) in Section 6, 
finding that the limits on and tan (3 are relaxed to 

> 46 GeV, tan/? > 1.9. (3) 

Section 7 is devoted to a discussion of the possibility of Higgsino dark matter in such a 
nUHM scenario. We find that the LEP II searches for charginos, neutralinos and Higgs 
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bosons together now exclude as dark matter an LSP that is more than about 70% Higgsino. 
We turn our attention to the Tevatron Colhder in Section 8. We compare the LEP II and 
Run I sensitivities to the MSSM parameters, and discuss and compare the regions of MSSM 
parameter space to which Tevatron Run II data should be sensitive [^. Finally, Section 9 
summarizes our conclusions and the prospects for future improvements and extensions of the 
analysis reported here. 



2 Theoretical Framework 

As already mentioned in the Introduction, we work in the context of the MSSM with R 
parity conserved. We assume a parameterization of soft supersymmetry breaking inspired 
by supergravity models with gravity mediation from a hidden sector. We assume that the 
soft supersymmetry-breaking scalar masses mo are universal at the supersymmetric GUT 
scale, as are the gaugino masses mi/2 and the trilinear parameters A. The renormalization 
of the physical values of the soft supersymmetry-breaking parameters is then calculated 



using standard renormalization-group equations. We use two loop RGEs|T^ to evolve the 



dimensionless couplings and the gaugino masses, and one loop RGEsp3| for the other soft 
masses, and we include one- loop SUSY corrections to /i|16[ and to the top and bottom 
masses [0 . 

Deviations from scalar-mass universality could easily be expected, for example in string- 
motivated models where their magnitudes could be controlled by flavour- dependent modular 
weights [0 . Upper limits on flavour-changing interactions place restrictions on the possible 
generation-dependences of scalar mass parameters, though these are relatively weak for the 
third generation. In any case, these do not constrain non-universalities between sparticle 
fields with different quantum numbers, namely vs £l vs vs ql- Nevertheless, we neglect 
such possibilities in our analysis. However, although our default option is that universality 
extends also to the soft supersymmetry-breaking contributions to the Higgs scalar masses 
(UHM), we do also allow for the possibility that their soft supersymmetry-breaking masses 
may be non- universal (nUHM). 

We use the renormalization-group equations and the one-loop effective potential to im- 
plement the constraints of a consistent electroweak vacuum parameterized by the ratio tan (3 
of Higgs vev's. We therefore adopt a parameterization of the MSSM in which mo,mi/2,^, 
tan (3 and the sign of are treated as independent parameters, with the magnitude of the 
Higgsino mixing parameter ^ and the pseudoscalar Higgs mass (or, equivalently, the bi- 
linear soft supersymmetry-breaking parameter B) treated as dependent parameters. In the 
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UHM limit, the correlation between mi/2 and /i is such that the LSP neutralino x typically 
is mainly a,U{l) gaugino B (Bino). Since /i becomes a free parameter (along with m^) in the 
nlJHM case, x "^^Y become Higgsino-like for certain parameter choices (roughly M2 > 2/i). 
Fig. |I| gives an overview of the /x, M2 = (a2/«G(7r) x ^1/2 plane for the illustrative choices 
tan/3 = 3, niQ = 100 GeV, At at its quasi-fixed point ~ 2.25M2 1113, and rriA = 1 TeV, 



showing various contours of the relic density the contour m^± = 100 GeV, contours 

of the mass of the lightest MSSM Higgs boson, and contours of Higgsino purity. In this 
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Figure 1: The fj,, M2 = {a2/ciGUT) x ^1/2 plane for tan/3 = 3, mo = 100 GeV and tua = 1 
TeV. Contours ofVL^h? = 0.025, 0.1 and 0.3 are shown as solid lines, and the preferred region 
with 0.1 < Q^h"^ < 0.3 is shown light-shaded. There are also dashed lines corresponding to 
m^± = 100 GeV. The near-horizontal dot-dashed lines are Higgs mass contours, and the 
hashed lines are 0.9 Higgsino and gaugino purity contours. The dark shaded region has 
m^± < mz/2. 

figure we have neglected neutralino-slepton coannihilation (discussed in detail below), since 
a small change in mo can move the masses out of the coannihilation region. We see that 
the LSP is mainly a i? in most of the fi, M2 plane displayed (where the relic density is of 
cosmological significance). One of the key questions we investigate is whether a Higgsino 
LSP is still allowed as a dark matter candidate by LEP II data P]. 

The relevance of direct LEP or Tevatron searches for sparticles does not need emphasis. 
In fact, it turns out that the indirect constraint on the MSSM parameter space provided 
by the Higgs search is also of great importance, as seen in Fig. |l|, particularly in the UHM 
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case. This constraint depends the MSSM mass parameters, because the mass of the hghtest 
MSSM Higgs h is sensitive, via radiative corrections, to sparticle masses, in particular the 
stop masses. This correlation has some impact even in the nUHM case, as we discuss in 
more detail later. 

It is interesting to confront the range of MSSM parameters still permitted by the LEP and 
other direct experimental searches for MSSM particles with other less direct experimental 
constraints, or with theoretical prejudices. Among the latter, one might mention gauge- 
coupling unification, lepton-quark mass unification and the absence of fine tuning. Although 
we consider all these prejudices appeahng, none of them is precise enough to enable us to 
draw any firm conclusions. Gauge-coupling unification cannot be used to constrain mo, mi/2 
and /i in the absence of a theory of GUT threshold effects. Lepton-quark mass unification is 
hostage to uncertainties in neutrino masses and mixing The fine-tuning price imposed 
by LEP data is rising [20|, particularly for small values of tan/5, but its interpretation is 
subjective and no consensus has been reached on the maximal pain that can be tolerated. 
The constraints we apply in this analysis are rather the indirect experimental ones provided 
by the measurement of 6 — > 57 decay whose implementation we discuss in the next 

Section, and the cosmological relic-density constraint already mentioned in the Introduction. 

In the parameter region of interest, the relic density ^^^^ increases with increasing 
^0i^i/2- Therefore the cosmological upper limit Vt^h"^ < 0.3 may be used here to set 
upper limits on these soft supersymmetry-breaking parameters. As discussed in the Intro- 
duction, strictly speaking there is no astrophysical lower limit on Q^h'^, even if one accepts 
that the cold dark matter (CDM) density QcDAih"^ > 0.1, since there might be other im- 
portant sources of CDM, such as axions or ultra-heavy relics. Nevertheless, one may take 
n^h^ > 0.1 ClS db default assumption. 

An important recent development has been the recognition that coannihilation of the 



LSP with next-to-lightest sparticles (NLSPs) may be important [|TOl in the Bino LSP region 
that is favoured in the UHM case, in particular. Generically, the NLSP in this region is the 
lighter stau fi, with the fin and not much heavier. Since the stable MSSM relic cannot 
possess electric charge, the allowed region of the MSSM parameter space is bounded by the 
line = irif^. Close to this line, x~ ^ ^"^^ I — I coannihilation effects suppress fi^/i^ below 
the range that would be calculated on the basis oi x~ X annihilation alone. This has the 
effect, in particular, of increasing the maximum allowable value of mi/2 and hence allowing 
< 600 GeV for f^^^^ — -^^^ larger mi/2 > 400 GeV, the allowed range of mo has 
a typical thickness 6mo ~ 30 GeV. On the other hand, when mi/2 ^ 400 GeV, there is a 
relatively broad allowed range for mo between about 50 and 150 GeV, depending on tan j3, A 
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and the sign of /i. 

We have shown previously |^ that the lower limit on imposed by data from LEP 
and elsewhere may be strengthened by combining it with additional theoretical constraints 
such as the cosmological relic density. The previous analysis included coannihilation effects 
only in the Higgsino region. The inclusion of LSP-NLSP coannihilation in the Bino region 
is less important for the inferred lower limit on m^, as we discuss later. However, it is 
important when one is assessing how much of the preferred range of MSSM parameter space 
may escape searches at LEP and elsewhere. 

In our analysis below, we consider parameter ranges that highlight the current experi- 
mental bounds and are consistent with the relic cosmological density. We use four default 
values for tan/5, namely 3, 5, 10, and 20. Lower values of tan (3 are disfavoured by the 
LEP Higgs mass limit, and the study of higher values would require an improved treatment 
of the cosmological relic density calculation for large tan (3, which lies beyond the scope of 
this paper. Although most of our Figures display more restricted ranges of mi/2, we note 
that in the UHM case cosmology allows values of mi/2 up to ~ 1400 GeV. We consider two 
possible treatments of the trilinear soft supersymmetry-breaking parameter Aq at the GUT 
scale which we assume to be universal. The conservative approach in the UHM case is to 
vary Aq so as to minimize the impact of the accelerator constraints (UHMmin), and the other 
is to choose [|I^] Aq = — mi/2 so as to maximize the area in the mo — rni/2 parameter plane 
in the present electroweak vacuum is stable, and CCB minima are irrelevant. In the nUHM 
case, we must also specify values of fi and m^. For the purpose of translating Higgs mass 
limits into limits in the itlq — mi/2 plane, we fix rriA = 10 TeV, so as to maximize the light 
Higgs scalar mass and therefore derive the most conservative bound possible. We allow fi 
and Ao to vary as much as possible while remaining consistent with the experimental lower 
bounds on the sparticle masses. 

There are restrictions on large values of Aq, so as to ensure that the sfermion masses are 
well-behaved. One of the most stringent bounds is that imposed by the experimental lower 
limit on the lighter stop mass, which depends on in the way depicted in Fig. |], which 



combines the constraints from p2| , p3| , p^ . Another important requirement is that the LSP 
not be a stau: m^^ > m^. The impacts of the stop and stau constraints are illustrated in 
Fig. Ij. We show for, /i > and tan/5 = 3,5,10 and 20, the corresponding upper limits 
on Aq in the UHM functions of mi/2 for mo = 100 GeV. Also shown in Fig. || as 

broken lines are representative Higgs mass contours. It is apparent that the Higgs mass is 
very sensitive to the value of Aq. The corresponding figures for ^ < are similar but allow 
somewhat higher values for Aq. The sensitivity of the Higgs mass to Aq translates into a 
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Stop - cx decays 




Figure 2: Present constraints in the {mi,m^) plane assuming a 100% branching ratio for the 
decay process ti— >m/cx- The vertical hatched hand represents the recent ALEPH exclusion 
valid for any AM value presented in jl^/; the light grey region is the LEP-comhined excluded 
region using data up to ^/s = 189 GeV under the most conservative assumption for the 
coupling Ztiti [2^]; the cross-hatched area is the CDF exclusion using the complete Run I 
data sample 
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Figure 3: Upper bounds on the trilinear supersymmetry-breaking parameter Aq, for tan/3 = 
3, 5, 10 and 20 and fi > , as a function ofmi/2. The shaded regions yield either a tachyonic 
i or a f LSP. Also shown is the dependence of the lightest Higgs mass on mi/2 and Aq. 
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3 Constraints from b ^ Decay 

The width for the inclusive decay B^Xs7 is determined by flavour- violating loop diagrams, 
and is therefore sensitive to physics beyond the Standard Model. In generic models with two 
Higgs doublets, significant contributions come from charged Higgs boson exchange, which 
always increases the SM prediction for = -B(B— >-Xs7), allowing severe lower limits on 
the mass of the charged Higgs boson to be set: see for instance. However, these limits 
do not apply directly to supersymmetric extensions of the Standard Model, because, in 
addition to the two Higgs doublets, there are chargino-stop contributions which can interfere 
destructively with the charged Higgs boson ones, and thereby reduce the predicted rate for 



Calculations including next-to-leading order (NLO) QCD corrections exist for both the 
Standard Model and generic two-Higgs-doublet models (see pSl and references therein). 



whereas in the case of supersymmetry the leading order (LO) calculations |27] have been 
complemented with NLO QCD corrections that are valid only under certain assumptions [28|. 

We include in our numerical analysis a Bgj calculation based on the full NLO treatment 
for the Standard Model and charged Higgs contributions, and the best available^ treatment 
of QCD corrections to the supersymmetric contributions The latter turned out to be 
of limited applicability in our analysis, since the conditions in which they well approximate 
the whole NLO supersymmetric corrections are usually not met. Therefore, the results pre- 
sented below are based mostly on the LO supersymmetric contributions only. The prediction 
for Bg-y depends on some experimental inputs and on three renormalization scales. The ex- 
perimental inputs are the top-quark mass, the CKM mixing-angle factor | VjfeVJ*/V"c6|, the c 
and b quark masses, the inclusive semileptonic branching fraction of B hadrons Bi^pt^Xi the 
strong coupling ^^(M^), and the electromagnetic coupling ctem- The renormalization scales 
are those relevant to the semileptonic and radiative processes ~ mb-, and the high-energy 
matching scale ~ Mw We used as nominal values and errors for these quantities those 
quoted in Table 1 of [^, except for Biept^x, Ois{Mz) and Oem- For the former we used the 



latest average provided by the LEP electroweak working group |2^, i?iept,x = 0.1058±0.0018 



For as{Mz) we took the latest PDG combination: 0.119 ±0.002. Finally, for aem we 
took the value at = following p8|. For each given point in the parameter space, we 



determined the theoretical error 5Bi^^"' as the RMS of 1000 S,^ values obtained by varying the 
experimental inputs with independent and Gaussian errors. Moreover, we determined the 



^The code implementing these calculations has been kindly provided to us by P. Gambino, who also 
helped in designing a recipe to determine the applicability of the supersymmetric NLO calculations. 
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reference theoretical prediction for Bl'^'^°^ conservatively, as the value closest to the measured 
one that we could obtain by varying independently the three renormalization scales from 
half to twice their nominal value. 

The experimental measurements of the rate for the inclusive process B— s>Xs7 s-re 
dominated by the latest CLEO result 

i?^^'^^ = (3.15±0.35±0.32±0.26) x 10-^ (4) 

which is in good agreement with the SM prediction of (3.29 ±0.33) x 10^^ |3l|]. There is 
therefore no need for any physics beyond the Standard Model, and we establish only upper 
limits on the possible supersymmetric contributions. 

To determine the 95% confidence-level exclusion domain in the supersymmetric parameter 
space, we treated separately the contributions to the CLEO error. The three terms in (3) 
come from limited statistics, experimental systematics and model dependence, respectively. 
By adding in quadrature the first two terms we defined 53^^°''^ which we treated as a Gaussian 

= n 9 

have then defined a function 



error, whilst we considered the third one, 63'^°'^'^'' = 0.26, as an additional scale error. We 



arsmeas ntheor I XfSmodel \ 2 



2 _ VI'-^S7 '-^57 I ^'-^sj J /gN 



and regard as excluded those points giving a probability for one degree of freedom smaller 
than 5%. 

We have investigated the impact of the Bs-y constraint (H, |^) in the (/i, M2) plane, without 
making the UHM assumption. Fig. ^ shows the domains excluded by Bs-y as a function of rriA 
for tan/5 = 3, mo = 500 GeV and At = 2mi/2, the quasi-fixed point value 0. The kinematic 
reach for charginos at ^/s = 204 GeV is also shown for comparison. As expected, the 



extent of the excluded region depends strongly on m^~ymj^± — m^, essentially vanishing 

when niA > 350 ^ 400 GeV. When = 350 GeV, the excluded region collapses to the 

cross-hatched area shown in Figs, ^a and Focusing on the case of = 250 GeV, at 

large M2 either the chargino or the stop is heavy enough to suppress the supersymmetric 

contributions to Bg-y- In this case, the positive charged-Higgs contribution dominates, making 

the predicted value of Bg-y incompatible with the measured value. For moderate M2 values, 

the supersymmetric contribution becomes sizeable. When /x > 0, it interferes negatively 

with the charged Higgs contribution, reducing the predicted value for Bg-y, whereas for /i < 

^Note that, in this Figure alone, we display plots for the renormalized low-energy value At, rather than 
the input GUT value Aq. 
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Figure 4: Constraints on the nlJHM parameter space imposed by Bg^: domains in the (/i, M2) 
plane excluded for tan P = 3 (a,b,c) and tan P = 10 (d). In all plots the 'reference' excluded 
region for rriA = 250 GeV, niQ = 500 GeV and the infra-red quasi-fixed-point value Aq = 
2mi/2 is shaded, assuming rrit = 175 GeV. The effect of varying niA is shown in panel (a), 
the effect of varying is shown in panel (b), the effect of changing the sign of A is shown 
in panel (c), and panel (d) illustrates the effect of increasing tan/?. Please see the text for 
further details. 
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it adds constructively to the charged Higgs contribution, strengthening the exclusion. This 
explains the shape of the excluded domains for niA = 150, 250 GeV. The excluded domains 
depend only mildly on mo, as shown in Fig. whilst it can be seen in Fig. ^ that the 
dependence on the sign of At is significant, and we also show for comparison the case At = 0. 
Finally, Fig. ^ shows the same exclusion domains as Fig. but for for tan/3 = 10. While 
the charged Higgs contributions essentially saturate for tan /5 > 4 -h 5, the supersymmetric 
contributions contain terms of order 1/ cos (3, and therefore increase with tan (3. When /i > 0, 
this has the effect of further reducing the prediction and hence the excluded region, while, 
for /i < 0, it enhances Bg^y, thereby extending the sensitivity of these constraints to larger 
rriA values. 



4 Update on constraints from LEP II 

The latest general presentations of results from the four LEP collaborations were made on 



Mar. 7th, 2000 [32]. They were based on the following mean integrated luminosities in each 



experiment at the indicated centre-of-mass energies: 



EcM 


= 188.6 GeV : 


171 pb" 


EcM 


= 191.6 GeV : 


28 pb" 


EcM 


= 195.6 GeV ■ 


78 pb" 


EcM 


= 199.6 GeV : 


80 pb" 


EcM 


= 201.6 GeV : 


38 pb" 



(6) 



No significant signals were announced in any sparticle or Higgs search channel. Numerical 
lower limits on the Higgs boson masses were presented separately by the four experiments. 



and a preliminary combination performed by the LEP Higgs Working Group is available [|T3 
Limits were also presented by the individual experiments on sparticle production within 
several frameworks, but the combination of the standard channels usually provided by the 
LEP Supersymmetry Working Group |l[] was not made available at this time. We extrapolate 
the available combined LEP limits, provided on the basis of the running up to Ecm = 
188.6 GeV for the sparticles and up to Ecm = 201.6 GeV for the Higgses, to include the 
higher-energy /luminosity data. 

We consider the following possible scenarios for the future evolution of the integrated 
LEP luminosity. The pessimistic one is that no significant additional high-energy luminosity 
is accumulated (remember the beer bottles?). In this case, the sparticle and Higgs sensitivity 
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will remain essentially as they are at the end of 1999. We believe that a more realistic scenario 
is for LEP to accumulate luminosity at the same average rate of 1.3 pb^^ as in 1999, but at 
somewhat higher energies, say 2/3 at Eqm = 202 GeV and 1/3 at Eqm = 204 GeV. This 
would result in the following total integrated luminosities per experiment at energies above 
EcM = 200 GeV: 

EcM = 202.0 GeV : 160 pb~^ (7) 
EcM = 204.0 GeV : 60 pb"^ 

A more optimistic scenario would be that luminosity is accumulated at a rate sometimes 
achieved in 1999, but not consistently, and that 50% of the running is at Ecm = 204 GeV 
and 206 GeV: 

Ecm = 202.0 GeV : 140 pb"^ 

Ecm = 204.0 GeV : 80 pb"^ (8) 
Ecm = 206.0 GeV : 20 pb"^ 

We do not provide detailed results for this optimistic scenario, but do make some comments 
on its potential impact. 

The sparticle final states of relevance for this analysis are X^X~ ^XX' yX'x' ^-^^d In 
addition, the experimental limits on the squark-production processes tt and bb can be used to 
infer constraints on the A parameters, as discussed in Section 3. We contrast two approaches 
in the following, either we take the CCB constraint into account and fix Aq so as to minimize 
its impact, or we take a conservative approach, allowing any value of A consistent with the 
experimental limits on m^^ and other constraints (UHMmin). 

The experimental efficiency for sparticle detection and hence the cross-section upper limit 
depends on other parameters besides the target sparticle mass, for example the mass differ- 
ence AM in the sparticle decay, e.g., x~^ ^ X + ^ ■ We have modeled these varying detection 
efficiencies using a multistep function, with a lower AM cutoff, low and high AM regions. 
We have used the available publications by the LEP collaborations and the documentation 
provided by the LEP Supersymmetry Working Group |1| to derive reasonable values for the 
transitional values of AM and the average efficiency values within the two regions. We did 
the same for the background contaminations, except in the case of slepton production, in 
which case we modeled the dominant W^W~ background in different regions of the (M^-, M^) 
plane using its detailed kinematics. We have checked that our parameterization reproduces 
the available published results. In each of the 'realistic' and 'optimistic' scenarios (|^, ^), we 
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LEP scenario 


1999 


'realistic' 2K 


'optimistic' 2K 


Max (GeV) 


201.6 


204 


206 




0.19 


0.11 


0.11 




0.06 


0.05 


0.06 




0.09 


0.08 


0.08 





~0.08 


~0.08 


-0.07 



Table 1: Examples of estimated upper limits on sparticle cross sections ( in pb ): for charginos 
and sleptons, assuming = 50 GeV, M^± = 100, M^± = 95 and Mf± = 90 GeV, respec- 
tively. In the case of neutralino production, typical values are given. 



make the assumption that the experimental efficiencies and contaminations remain similar 
to those at Ecm < 189 GeV, based on the fact that the properties of the standard process 
do not change dramatically in the spanned energy range. Examples of the estimated upper 
limits on sparticle production cross sections that we obtain from our extrapolation to the 
three LEP running scenarios discussed above are given in Table 1. 

For charginos, we conservatively assume no detection efficiency for AM < 5 GeV. For 
larger AM values, the estimated upper limits allow one to exclude chargino production up 
to a few hundred MeV below the kinematic limit, unless sneutrino masses, and hence mg, are 
very small. In the case of associated neutralino production, we combine all the kinematically 
accessible channels, weighted by their visible cross sections, i.e., we take into account the 
branching fractions into final states, and the estimated efficiencies. We found that the 
estimated upper limits depend only weakly on the point in the MSSM parameter space, and 
typical values are given in Table 1. In the case of slepton production, we set the AM cutoff 
at 3 GeV. As an example. Fig. |^ shows the exclusions we obtain in the plane {mg^,m^) for 
tan/5 = 3, /i = 200 GeV and BR(e^ — > xe^=100%), under the different hypotheses for LEP 
running in 2000. 

We now turn to the estimation of upper limits on Higgs production. In order to estimate 
the prospects for Standard Model Higgs limits from the reaction e~^e^ — » hZ^, we take the 
simple parameterization of the LEP limits obtained from data up to Ecm ~ 189 discussed 
in [^. Extrapolating this parameterization to include the 1999 data set (|^) leads to the 
estimated limit mn > 109 GeV. This result is in good agreement with the expected limit 



reported in [|13l. However, the observed limit quoted in the same reference is 107.9 GeV: 
the difference is explained as a statistical fluctuation in the data at the level of one standard 
deviation. 



15 



"4 - |i= 200 GeV, tan|3= 3 




95% C.L. exclusion 



^ — ^ — , — ^ — ^ — ^ — ^ — , — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ — ^ 

50 60 70 80 90 100 

m„ (GeV/c ) 

eR 



Figure 5: Constraints in the mg^,m^ plane imposed by the combined LEP data at Ecm < 
189 GeV (dotted line), our estimates for the limits obtainable by combining the LEP data 
taken in 1999 ^ (dashed line), and our estimates for the possible 2000 exclusions in the 
'realistic' (^ (solid line) and 'optimistic' (dot-dashed line) scenarios described in the 
text. 



The corresponding limits for the 'reahstic' and 'optimistic' running scenarios for the 
year 2000 are 112 and 114 GeV, respectively, following roughly the empirical rule Mh > 
-Ecm— 93 GeV. Similar estimates apply to the MSSM for small tan /? < 5, as shown in Fig. 
For larger values of tan/5, we use the same limiting cross section for e'^e^ hZ^, which 
gives a weaker lower limit on m^, because of the smaller Z^Z^h coupling. When tan/? > 
8, the production mechanism e'^e" HA becomes important, which we include in our 



analysis following again the prescription given in ||33[. Our estimated limiting curves in the 

m/j,tan/3 plane shown in Fig. ^have been calculated in the 'Max(M/j)' benchmark scenario 

suggested in . We have not attempted to combine the hZ^ and HA analyses, but have only 

overlapped them, so our results could be considered conservative in the intermediate-tan /3 

■^We have verified that these hmits are not weakened by the appearance of invisible decay modes h, A ^ 
XX- 
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region. We indicate in Fig. |^ the LEP limit for the full 1999 data set (dot-dashed line), our 
estimate for the 'realistic' 2000 running scenario (|^) (shaded) and the 'optimistic' scenario 
(|) (dashed line). Also shown (shaded) are the regions of the rrih, tan P plane excluded by 
theoretical calculations [03 . 




Figure 6: Constraints in the ruh, ta.n (3 plane imposed in the 'Max(Mh)' benchmark scenario 
by combining the LEP data taken in 1999 ^ (dot-dashed line), and our estimates for the 
possible 2000 exclusions in the 'realistic' ^ (solid line) and 'optimistic' ^ (dashed line) 
scenarios described in the text. Also shown (dark-shaded) are the regions of the nih, tan P 
plane excluded by theoretical calculations 

Some caution is required when we compare our results directly with the lower limits given 
by the LEP experiments because of theoretical uncertainties in the MSSM Higgs mass cal- 
culations. Conservatively, we allow for an error of ~3 GeV in these, so that we translate the 
experimental limits into the supersymmetric parameter space using the theoretical contours 
for 104 GeV (for tan f3 = 3) and 100 GeV (for tan (3 = 5) in the case of the complete 1999 
data LEP scenario, and 109 GeV and 108 GeV, respectively for the two values of tan/3, in 
the 'realistic' 2000 LEP scenario. At the higher values of tan (3 considered, the Higgs mass 
limits do not provide strong constraints and are not used. 
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We stress that radiative corrections to chargino and neutralino masses, though less dra- 
matic than to Higgs masses, are also relevant to the interpretation of experimental limits on 
physical particle masses in terms of constraints on MSSM parameters such as /i, mi/2 and 
mo 0. Two such effects are seen in Fig. 0. We see that the differences between the 1999 and 
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Figure 7: The effects of radiative corrections to chargino masses in (a) the (/i, M2) plane 
for tan (3 = 5 and (b) the {mi/2,mQ) plane for tan f3 = 5 and fi> 0. The contours obtainable 
from the 1999 and 'realistic' 2K data are indicated, both with (thicker lines) and without 
radiative corrections (thinner lines). We notice that the differences between the lines with 
and without radiative corrections are larger than those between the 1999 and 2K lines. 

2K limits are considerably smaller than the shifts induced by the radiative corrections. In 
particular, as shown in the /i, mi/2 plane in panel (a), radiative corrections are very signifi- 
cant in the delicate Higgsino region discussed in Section 7. Their inclusion is indispensable 
for the accurate interpretation of the LEP data, as we do throughout this paper. 

5 The Case of Universal Higgs Masses 

We next apply the above accelerator constraints under the assumption that the soft super- 
symmetry breaking masses are universal, including the Higgs multiplets (UHM), exploring 
their impact in the the mi/2,'"^o parameter plane and comparing them with the constraints 
from cosmology on the relic abundance of the LSP. We remind the reader that, in this UHM 
context, for fixed tan/3 and sign of /i, the only parameter choice remaining is the value of 
Aq. We discuss below two cases, one in which we require the absence of charge and color 
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breaking (CCB) minima |]12|, but fix Aq = — mi/2 so as to minimize their impact, and the 
other in which we disregard CCB minima, and allow Aq to vary freely (UHMmin). 

We start with the CCB UHM case shown in Figs. ^ and ^. One is safe from CCB minima 
above the curved solid lines in these plots, which are calculated with Aq = —mi/2, so as to 
minimize the impact of this prospective constraint. Plotted as near- vertical dashed lines are 
the chargino mass contours: = 102. At the higher values of tan/3, the bound on mi/2 
from the chargino mass limit is nearly independent of tan (3, and the dependence is always 
very slight for fi > 0. For completeness we also show the limit from the selectron mass 
bound. 

The light shaded regions in Fig ^ are those excluded by the b s'j constraint discussed 
in Section 3. We see that, for fi < 0, the impact of Bg^ constraints increases sizeably with 
tan /3, and covers in these cases a significant fraction of the region otherwise preferred for 
dark matter reasons. The cut-off at large mo, mi/2 is due to the corresponding increase 
of rriA, and hence mH±, in the UHM, which reduces the charged Higgs contribution. The 
exclusion is not very sensitive to the Aq value chosen. As we can see in Fig |^, for /i > the 
interference between the supersymmetric and charged Higgs contributions cancel the effect 
of new physics in the low-medium tan (3 range; however, there is still some sensitivity at large 
tan j3 where the large negative supersymmetric contribution make significantly smaller 
than the measured value, as seen in panel (d) of Figs. 

Also shown in these figures by near-vertical dot-dashed lines are the limits coming from 
the Higgs mass bounds. Note that for yU < 0, these contours only appear for tan/9 = 5, 
where we display the 100 and 108 GeV contours corresponding to the 1999 and prospective 
'realistic' 2K experimental limit, allowing a safety margin of 3 GeV as discussed earlier. At 
tan/5 = 3, the position of these contours is far off to the right, excluding the entire region 
displayed. At tan j3 = 10 and 20, the contours would appear to the left of the chargino 
bound and are not shown. For /i > 0, the limits are weaker, i.e., the contours move to the 
left. In the case of tan/5 = 3, the 102 GeV and 104 GeV contours have now moved into the 
displayed range of mi/2, and the contours for 100 GeV and 108 GeV are shown for tan (3 = 5. 
In Fig. [ly, we show an extended range in mi/2 and the position of the 104 GeV contour for 
tan P = 3, fi > 0. 

The Higgs mass contours depend on mg, because the radiative corrections to the Higgs 
mass cause these curves to bend left at very large rriQ. Thus ultimately, the only pure 
accelerator constraint on mi/2 comes from the chargino mass limit. However, as we have 
demonstrated previously, cosmology excludes such high values of mo, thus maintaining the 
importance of the Higgs mass bound in limiting mi/2, and ultimately m^ and tan/3. The 
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Figure 8: The mi/2, mo plane for < 0, A = — mi/2 so as to minimize the impact of the CCB 
constraint (indicated by a solid line) and tan/? = (a) 3, (h) 5, (c) 10 and (d) 20. The region 
excluded by our b —>■ s'j analysis has light shading. The region allowed by the cosmological 
constraint 0.1 < ^x^"^ — O-^; after including coannihilations , has medium shading. Dotted 
lines delineate the announced LEP constraint on the e mass and the disallowed region where 
ruf^ < has dark shading. The contour m^± = 102 GeV is shown as a near-vertical dashed 
line in each panel. Also shown as dot-dashed lines are relevant Higgs mass contours. 
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Figure 9: The plane for /i > 0, A = — mi/2 and tan/3 = (a) 3, (h) 5, (c) 10 and 

(d) 20. The significances of the curves and shadings are the same as in Fig. The light- 
shaded region in panel (d) is excluded by the b s'j constraint. The long dashed curves in 
panels (a), (b) and (c) represent the anticipated limits from trilepton searches at Run II of 
the Tevatron j^. 



21 



medium-shaded regions in Figs. ^ and ^ show the areas in the mo,mi/2 plane for which 
the rehc cosmological density falls between 0.1 < Qh"^ < 0.3 when co-annihilation effects are 
included. We note that the chargino mass constraint now essentially excludes the re-entrant 
parts of the dark matter density contours caused by resonant direct-channel annihilations 
when mi/2 < 160 GeV, which were visible in Fig. |I| as well as Fig. ^. The dark shaded 
regions in Figs. || and || correspond to a charged LSP, as indicated. 



We show in Fig. ^ the 'tail' of the cosmological region where 



for tan/5 = 3. 



As can be seen in Fig. 7 of [|T3], the tip of this region is allowed by the CCB constraint for 
tan /3 = 3, 10, and we can see in Fig. ^ that the CCB constraint is weaker for tan [3 = 5. We 
show in Fig. 10 the rrih = 104 GeV contour, corresponding to the 1999 bound on the Higgs 
mass after allowing for a 3 GeV theoretical uncertainty in the prediction. We recall that the 
Higgs mass limit is even stronger for negative fi. We can safely set the limit tan/3 > 2.8 for 
yU > in the UHM. Overall, the limits we obtain on tan/5 in different LEP scenarios for both 
signs of /i are shown in Table 2. 





1999 


'realistic' 2K 


/i < 
/i > 


3.2 
2.8 


4.0 
3.6 



Table 2: Limits on tan/3 imposed in the UHM by the 1999 and 'realistic' expected 2K Higgs 
mass limits. 



We now repeat the above UHM analysis for the more conservative case (UHMmin) in 
which we do not require stability against collapse into a CCB vacuum, and we allow Aq to 
vary as far as possible, consistent with the experimental constraints on m^^ in particular. 
Fig. |ll] shows options for tan (3 and the sign of which exhibit interesting differences from 
the previous CCB case. For simplicity, we have chosen not to include the regions that are 
excluded by the 6 — > S7 constraint: these turn out to be essentially independent of Aq, and 
hence may be taken from the corresponding panels in Fig ^ The interesting and significant 
differences are in the contours of the Higgs mass. The Higgs mass is sensitive to Aq, and 
may be significantly lower than in the previous CCB case, with corresponding implications 
for the lower limits on and tan (3 that we quote below. 

Clearly, by allowing Aq to vary (rather than restrict its value to —"^1/2), we expect weaker 
bounds from the Higgs mass than those found when the CCB constraints were incorporated. 
Since the Higgs mass constraint was only important for lower values of tan /3, we show results 



only for tan/3 = 3, 5 in Fig. y. Whilst, for tan/3 = 3 and < 0, the 104 and 108 GeV 
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Figure 10: An extension of the mi/2,'"^o plane for n > 0, A = — mi/2 and tan/3 = 3. Below 
the solid diagonal line, the LSP is charged and hence excluded. The absolute upper bound to 
mi/2 is found when the shaded region drops entirely below this contour. 

Higgs mass contours are still to the right of the displayed region in the figure, we see that 
in the other cases shown, all of the contours are moved substantially to the left. In fact 
for tan/? = 5 and /i > 0, the Higgs mass bound is no longer competitive with the chargino 
bound. As in the previous UHM case, we also find lower bounds on tan/5 in this case where 
the CCB constraint is relaxed, as shown in Table 3 for different LEP running scenarios and 
the two signs of fi. 





1999 


'realistic' 2K 


/i < 
/i > 


2.7 
2.2 


3.1 
2.7 



Table 3: Limits on tan/3 imposed in the UHM by the 1999 and 'realistic' expected 2K Higgs 
mass limits, relaxing the requirement that there be no CCB vacuum. 



6 Bounds for Non- Universal Scalar Masses 

In the previous Section, we have derived stringent limits on the mi/2,mo plane from the 
absence of sparticles and Higgs bosons at LEP, assuming universality for scalar masses 
including the soft Higgs masses (UHM). These limits are particularly strong at low values 
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Figure 11: The mi/2,mo plane for /i < and (a) tan/5 = 3, (h) tan/3 = 5, and /i > 
and (c) tan/5 = 3, (d) tan (3 = 5, with Aq allowed to vary, and no CCB constraint 
applied. For clarity, the region excluded by the b —* S'-y constraint has not been shaded: it 
is essentially identical to that in Fig. |^. The region allowed after including coannihilations, 
by the cosmological constraint 0.1 < ^x^"^ — ^-^ has medium shading. Dotted lines delineate 
the announced LEP constraint on the e mass and the disallowed region where mf^ < has 
dark shading. The contour m^± = 102 GeV is shown as a near-vertical dashed line in each 
panel. Also shown as dot-dashed lines are relevant Higgs mass contours. 
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of tan/3, and in fact exclude tan (3 < 2.8. One should expect these limits to weaken when 
the assumption of UHM is relaxed (nlJHM). In this Section, we rederive the appropriate 
limits in the mi/2, ttlo plane for the more general nUHM case. We now treat both fi and 
as independent parameters, in addition to the free parameters mi/2,mQ, A, and tan/5 of the 
previous Section. As indicated in Section 2, we take rriA = 10 TeV, so that the Higgs mass 
limits give the most conservative bounds on mi/2- As before, we restrict the values of A by 
requiring that m^^ be consistent with the experimental lower limit, and > m^^, as shown 
in Fig. H. 



Our results for the nUHM case are shown in Figs. |12| and |T3|. We again show the 
kinematical limit on the chargino mass: = 102 by the near-vertical dashed line. Again, 
for /i > 0, these lines are essentially vertical, because the chargino mass is independent of 
mo, apart from the effects of radiative corrections. As seen in the different panels of Fig. |12|, 
the lower bound on mi/2 from the chargino bound increases from mi/2 = 112 GeV to 145 
GeV as tan/5 is increased from 3 to 20 for fi < 0. In contrast, as seen in Fig. 13, the bound 
on mi/2 for /i > lies near 140 GeV over the same range in tan/?. As expected, the nUHM 
curves for the chargino mass limits always lie to the left of the UHM ones. 



As in Figs. ^, ^, and 11, the shaded region corresponds to the parameter values in which 
it is possible to achieve 0.1 < Qh'^ < 0.3. This region is unbounded from above, since it 
is possible to adjust fi to insure an acceptable relic density even if the sfermion masses are 
large. This is because by lowering /i, the LSP can become a mixed state (rather than an 
almost pure bino) and annihilation channels via Z°-exchange open up. 

We also show in Figs. [l2| a (b) and |l^ a (b) the nlJHM contours for Higgs masses of 104 
(100) and 109 (108) GeV, the 1999 and 2K 'realistic' bounds for tan/5 = 3 (tan/5 = 5). For 
tan/3 = 5 and yU > 0, the 100 GeV contour is to the left of the chargino mass contour and 
is not shown. We recall that, in the nlJHM case, there are no unique Higgs mass contours 
in the mo,mi/2 plane, due to the freedom in choosing /i, m^i and Aq. The contours shown 
correspond to parameter choices giving the weakest bound. 

Because the Higgs mass bound is weakest for relatively large values of fi in the nUHM 
case (implying that the neutralino is a bino), the relic LSP density increases with ttlq as one 
moves upward along the Higgs mass contours and the sfermion masses increase. In some 



cases, e.g., that in Fig. |T3|a, the relic density along the 104 GeV contour exceeds 0.3 at 
about mo = 140 GeV. At higher values of mo, the value of mi/2 must be increased to remain 
consistent with both the Higgs mass limits and cosmology. This adjustment is shown by the 
dotted curve to the right of this contour. Similar behavior was seen in |21|]. For the other 
Higgs mass contours, either the shift (in mi/2) is insignificant at mo < 200 GeV, or the relic 
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Figure 12: The mi/2,mo plane for fi < in the nUHM case, for tan/3 = (a) 3, (b) 5, 
(c) 10 and (d) 20. The region allowed after including coannihilations, by the cosmological 
constraint 0.1 < Vt^h"^ < 0.3 is shown shaded. Dotted lines delineate the announced LEP 
constraint on the e mass and the disallowed region where rrir^ < has dark shading. The 
contours m^± = 102 GeV are shown as near-vertical dashed lines. Also shown as near- 
vertical dot-dashed lines are Higgs mass contours. 
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Figure 13: The mi/2, mo plane for /i > m the nUHM case, for tan/3 = (a) 3, (h) 5, (c) 10 
and (d) 20. The significances of the curves and shadings are the same as in Fig. |7^. The 
dotted extension of the mu = 104 GeV contour corresponds to the shift in the Higgs contour 
when the cosmological limit on the relic density is imposed (visible only in (a)). 



27 



density does not exceed flh"^ = 0.3 for niQ < 200 GeV. For /i < 0, the relic density is never 
saturated for the values of mo shown. As before, for tan/3 = 10 and 20, the chargino bound 
is always stronger than the Higgs mass bound. 

7 Higgsino Dark Matter 

We now turn to the question whether there is any room left for Higgsino dark matter. We 
update the analysis of P], including the improved experimental limits discussed above, and 
we explore the sensitivity of our conclusions to the possible range of LSP relic density. We 
begin by reviewing briefly the analysis of 0. 

We recall that a general neutralino is a linear combination of the Higgsinos and neutral 
gauginos, x = + aW^ + •yHi + 6H2. In this notation, the Higgsino purity is defined to be 
p = Vo^~+~^! 3' state that is half gaugino and half Higgsino has Higgsino purity 1 / -\/2, and, as 
in 0, we take as our working definition that a neutralino is a 'pure' Higgsino if p > 0.9, even 
though such a state already has a sizeable gaugino fraction: y^a'^ + (3'^ ~ 0.44. The lightest 
neutralino tends to be Higgsino- like if fi < M2/2 and gaugino-like when fi > M2/2. This was 
already shown in Fig. |1|, where we plot as hashed lines contours of Higgsino purity in the 
/i,mi/2 plane for tan/3 = 3. We have also plotted thin solid contours for Q^h"^ = 0.025,0.1 
and 0.3, thick solid lines corresponding to m^± = 100 GeV, and dot-dashed Higgs mass 
contours. In this illustration we have taken rriA = 1 TeV and mo = 100 GeV. It is apparent 
that the bulk of the cosmological region with 0.1 < ^x^"^ — ^■'^ larger \^\ (for given 
M2) than do the Higgsino purity contours, indicating that LSP dark matter is generically a 
gaugino: in these regions, it is mainly a Bino. There are, however, small regions at smaller 

(for given M2), where the LSP is mainly a Higgsino. However, as can be seen in Fig. |l|, 
this Higgsino possibility is under severe pressure from several LEP constraints, including the 
chargino, xx' ^^^^d Higgs searches, and it is the possible exclusion of the Higgsino dark matter 
regions we explore in this Section. 

Accordingly, we now focus in more detail on the Higgsino regions, as illustrated in Figs. 113 



and 15, where detailed views of the Higgsino parts of the fx, M2 plane are shown for tan/3 



2, 3, 5 and 10. Consider in particular Fig. 0a, for tan/5 = 2 and /i < 0. Here we have taken 
TTiA = 10 TeV to minimize the effect of the Higgs mass limit, and mo = 1 TeV to maximize 
the neutralino relic density. In contrast to Fig. |1|, we have also adjusted At to maximize the 
Higgs mass and produce the weakest Higgs constraint. The hashed, dot-dashed, thin solid 
and dark thick solid contours are as in Fig. We plot as a dashed line the current chargino 
mass limit. We also show as two solid contours the most recent 1999 LEP 2 bounds on the 
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summed visible cross section for associated neutralino production cr(e+e^ XiXj)vis and 
our 'realistic' estimate for the final 2K bounds (see Table |I|). We recall that the associated 
production bounds are more constraining for smaller values of the scalar masses 0. It is 
evident that the entire region with ^^^^ > 0.1 and Higgsino purity p > 0.9 is excluded by 
the current experimental limits, for this value of tan (3. 

The general shape of the f^^^^ contours can be understood as follows. As one moves to 
large values of M2, the neutralino becomes more pure Higgsino, which leads to an approxi- 
mate three-way mass degeneracy between the lightest and next-to-lightest neutralinos X2 
and the lightest chargino x^- Since the X2 and x^ are therefore abundant at the time when 
the X freezes out of chemical equilibrium in the early Universe |^|, their coannihilations with 
the X act to bind the x more tightly in chemical equilibrium with the thermal bath, and 
delay the freeze-out of the x relic density. Since the X2 and x^ annihilate very efficiently. 



this greatly reduces the relic density of neutralinos for larger values of M2 [36|. We have 
included the one-loop radiative corrections to the chargino and neutralino masses, which can 
significantly affect both the experimental chargino limits and the neutralino relic density, 
when the chargino and neutralino are closely degenerate. Whenever the mass degeneracy is 
sufficiently tight to affect the chargino bounds, coannihilation suppresses the relic density to 
very small values, below those of cosmological interest [0]. 

Similarly, as one moves to larger values of the mass of the x increases, until > 
mw- At this point, the x can annihilate efficiently into W pairs, and the relic density 
drops dramatically as one crosses this threshold [|. The light solid contour corresponds to 
= m\Y, and the falloff of the relic density above this threshold is evident. It is the 
dramatic decrease in ^^^"^ ^^ir > my/ that not only excludes Higgsino dark matter, but 
also implies that we are not very dependent on our default choice of lower relic density cutoff: 
Vtyji^ > 0.1. In fact, one cannot even supply enough neutralinos to provide the galactic 
dark matter: Q^^h'^ ~ 0.025 above the W threshold, while still satisfying the experimental 
constraints. We see the same effect in all four panels of Figs |14|, and, although we only 
display four values of tan (3, we have verified that this is true for all tan (3. 

The situation is similar for /i > 0, as seen in Fig. In this case, the dominant experi- 
mental constraint comes from the chargino limits, which are shown as dashed lines for both 
1999 and the 'realistic' 2K scenario, whereas we plot only the 1999 contour for e'^e^ —* XiXj- 
The chargino constraints alone exclude f^y/i^ > 0.025 for a Higgsino-like neutralino. In all 



''For the sake of exposition, we forget for the moment that, at this low value of tan (3, the entire displayed 
region has a Higgs mass less than 106 GeV and can be excluded on this basis alone. 

^Sub-threshold annihilation into W pairs smoothes out this sudden drop |35 , and shifts the left edge of 



the n^h^ — 0.1 contour a few GeV to the right 
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Figure 14: Small regions of the /i, M2 plane for tan/3 = 2,3,5,10 and = 10 TeV. 
The thin solid lines are contours for ^^^"^ — 0.1,0.5,0.025, and the light-shaded region is 
cosmologically preferred. The dashed (light solid) lines correspond to m^± = 100.3 GeV 
(m^ = mw ), dot-dashed horizontal lines correspond to the indicated Higgs masses, and the 
two darker near-vertical solid contours indicate the current neutralino associated production 
bound a^e^e" XX2,---)vis o-nd our estimate of the 'realistic' final hound (see Table |7p. 
Hashed contours represent Higgsino purity. In panel (a) [(d)], the Higgs mass is everywhere 
less than 106 GeV [greater than 109 GeV]. The dark-shaded regions in panels (c) and (d) 
are those surviving all the constraints. 
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cases, an interesting amount of cold dark matter is possible only if < 0.7, i.e., the LSP is 
either predominantly a gaugino or a strongly mixed state. We conclude that a predominantly 
Higgsino state cannot provide a substantial component of the dark matter. 

This conclusion is robust with respect to variations in the other sparticle masses. The 
sfermion masses have already been taken large enough for the contribution to neutralino anni- 
hilation from sfermion exchange to be negligible, and, as already noted, lighter sfermions yield 
tighter constraints from associated neutralino production. The experimental chargino limits 
fall below the kinematic limit when the sneutrino is closely degenerate with the chargino: 
however, again, a light sneutrino enhances the associated neutralino production limits, so 
this also provides no loophole. A cosmologically interesting relic density can therefore only 
be achieved by either heavily mixed or pure gaugino neutralino states. 

8 Limits from the Tevatron Collider 

The limits on the mSUGRA (UHM) parameter space from Run I of the Tevatron are sum- 
marized in 0, and we have already made use of their lower limits on stop masses. Two other 
DO limits are reported in 0: one is for squark and gluino jets and missing energy, and the 
other is for dilepton events. The former analysis is for tan/? = 2, and hence is not directly 
comparable with our plots. The analysis also assumes Aq = and /i < 0, and yields a lower 
limit 

mi/2 > 100 to 50 GeV (9) 

for mo between and 300 GeV. The dilepton analysis in has been performed for several 
values of tan/5 including the values 3 and 5 studied in this paper. Again for Aq = and 
/X < 0, the lower limit 

mi/2 > 60 to 40 GeV (10) 

was found for mo between and 300 GeV. These Tevatron Run I limits do not constrain 
the MSSM parameter space as strongly as do the chargino and Higgs limits from LEP II, 
unless scalar mass universality between squarks and sleptons is relaxed. We have therefore 
not displayed the Run I bounds in Figs. | and 

On the other hand. Run II of the Tevatron, scheduled to begin in 2001, will impose strict 
new limits at low my2- The dominant experimental constraint is expected to come from 
searches for trilepton signatures of xfx2 production |^. We display in Fig. ^ as long- 
dashed contours the anticipated reach of the upgraded Tevatron in this channel, assuming 
2 fb~^ integrated luminosity. We have taken the published curves from p|, which does not 
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display results for tan/5 = 20 or < The bounds are tightest at low mo, where they 
cut into the cosmological region as shown. At tan/3 > 3, the Run II curves bend over and 
intersect the LEP2 chargino bound within the cosmological region. Their most significant 
impact is for low values of tan/3. However, since the LEP II Higgs bounds dominate at 
low tan/3, the trilepton analysis at the Tevatron Run II will therefore not increase the 
absolute lower bound on the neutralino mass given by our analysis in this region. However, 
if the Tevatron is able to improve the LEP Higgs bound, this could raise substantially the 
neutralino limit at low tan/3. The Higgs bound is no longer important for tan/3 = 10, but 
in this case the Tevatron limit bites away less of the region favoured by cosmology. 



9 Summary and Conclusions 

One of our principal goals in this paper has been to obtain strengthened lower limits on 
the neutralino mass, combining the latest LEP data with the cosmological dark matter 



requirement 0.1 < il^h < 0.3. We summarize our limits in Figs. |T6|, under various differ- 
ent assumptions: universal (UHM) or non-universal (nUHM) soft supersymmetry-breaking 
scalar masses for Higgs bosons and (in the former case) whether one requires the present vac- 
uum to be stable against transition to a charge- and colour-breaking (CCB) vacuum or not 
(UHMjnin). Also, we give limits both for the available 1999 LEP data and with a 'realistic' 
assessment of the likely sensitivity of data to be taken in 2K. 

In all cases, for both positive and negative /x, the lower limits on are relatively 
insensitive to tan /3 at large tan j3. Here they are determined by the LEP chargino bound, 
as the LEP Higgs mass bound is weaker than the chargino bound at large tan/3. In fact, 
in the two UHM cases shown, the points at which the limiting curves bend upward, as 
one decreases tan j3, are precisely the points at which the Higgs mass bound becomes more 
stringent than the chargino bound. In the UHM cases, the neutralino mass limits are strong 
at intermediate values of tan /3 ~ 4-7 because, as discussed earlier, the cosmological bound 
on the relic density prohibits going to large values of mo, and ensures that the Higgs bound 
places a strong constraint. Below this break point, the lower limit on m^ increases rapidly 
with decreasing tan j3. Above this break point, the limit on m^ is relatively insensitive to the 
additional theoretical assumptions made, such as UHM vs. UHMmin or nUHM. However, in 
the nUHM cases, because one can increase mo sufficiently to weaken the Higgs mass bound, 
the break point occurs at a lower value of tan j3. To go to lower values of tan /3 then requires 



^The tan/3 = 3, 10 Run II curves and the tan/3 = 5 Run II curve in Fig. ^ come from separate analyses 
and reflect slightly different confidence levels, 99% and 3(t, respectively. 
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a substantial increase in m^. 



n 200 

1999 nUHM 
1999 UHM 

1999 UHM^i„ " 

2000 nllHM 
2000 UHM " 
2000 UHM^i„ ' 

- c'^ioo 

- - ' 50 



a) jjL< 



5 6 7 8 9 10 

tarijS 



20 



■ 


1 > II 


III 1 

1999 nUHiM 






1999 UHM 


I '■ I '■ 


\ 
\ 


1999 UHM^.^ " 




1 

. ; \ 


2000 nlJHM 




\ 

■ ^ 


2000 UHM 


1 '■■ \ 
1 '■■ \ 

: W 

1 ■, 
1 

1 

1 

1 

L— 


■ 1 \ 
': \ 1 

\ '. \ \ 
\ ■ \ V 
\ ■■. \ \ 

■■. \ \ 

\ \ \ 
■■■ \ \ 


2000 UHM^,^ " 


1 


1 1 1 


b) l^l> . 

Ill 1 



6 7 8 9 10 

tanjS 



20 



Figure 16: Lower limits on the neutralino mass as functions o/tan/3 for (a) fi < and 
(h) /X > 0. The curves correspond to the final 1999 LEP results (thin lines) and our 'realistic' 
expectations for the 2K LEP run (thick lines). We show the UHM case with Aq = —mi/2 
to avoid CCB minima (dashed curves): these are the strongest constraints. We also show 
(dotted lines) the more general UHM^i^i case where Aq is left free, and we do not require 
the absence of CCB vacua. We also display additionally the nUHM case, which is the most 
conservative and allows both ii and tua to be free in addition to Aq . 

In the UHM cases with and without the restriction forbidding CCB vacua, the lower 
limit on the lightest MSSM Higgs mass, in particular, implies lower limits on tan/5 which 
are plotted in Fig. The limits are somewhat stricter for fi < than for /i > 0, whether 
(UHM) or not (UHMmin) one requires the absence of CCB vacua. Indeed, Fig. Owill enable 



the appropriate conclusion to be drawn from whatever lower limit on the Higgs mass LEP 
eventually provides. We recall that the existing Higgs mass calculations in the MSSM are 
believed to be accurate to about 3 GeV. Therefore in computing the bounds on tan (3 for 
Tables |^ and |^, for example, we have conservatively shifted the exclusion curves of Fig. ^ by 
3 GeV to the left before reading the values of tan (3 off of Fig. We also show in Fig. [l^ 
the lower bound on tan f3 obtained in the nUHM, which is significantly weaker than in the 
UHM cases, and essentially independent of the sign of fi. 

If LEP does achieve the 'optimistic' 2K energies and luminosities (P), the above con- 
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Figure 17: Lower limit on tan/? imposed by the experimental and cosmological constraints, 
as a function of the experimental Higgs mass limit. The UHM, UHMmin o,nd nllHM labels 
are as in Fig. The > curve in the nUHM case is very similar to the fi < curve. 

straints will be somewhat tighter. The horizontal segments of Fig. [1^, corresponding to the 
chargino limits, increase by a fraction of a GeV; the vertical branches move to the right, 
intersecting the horizontal segments at tan/3 = 8(7.5) for < (/x > 0). And lastly, the 
lower limits on tan /3 improve to 





UHM 




nUHM 




4.7 


3.4 


2.3 


fi>0 


4.2 


3.0 


2.3 



Table 4: Limits on tan f], assuming the 'optimistic ' 2K energies and luminosities 

In many respects, LEP has provided the most stringent constraints on the parameters of 
the MSSM. This is true, in particular, for its lower limits on the Higgs mass, and the chargino, 
neutralino and slepton constraints from LEP compare favourably with the Tevatron bounds 
on squark and gluino masses, once the different mass renormalizations of electroweakly- 
and strongly- interacting sparticles are taken into account. As we have shown, the present 
LEP data may be combined to set interesting lower bounds on the lightest neutralino mass. 
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in particular if it is assumed to constitute the dark matter favoured by astrophysics and 
cosmology, and on tan (3. It may well be that these lower limits will be further strengthened 
by the LEP run during 2000, as we have discussed in this paper. However, this would be 
the pessimistic scenario. There is still a chance that sparticles or the Higgs boson may turn 
up this year, in which case we would be delighted to see our bounds superseded. LEP may 
not yet have discovered supersymmetry, but it certainly deserves to! 
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